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Queues arising in buffers due to either random interruptions of the 
channel or variable source rates are analyzed in the framework of a single 
switched system. Examples of systems to which the results of the paper 
may be applied are: multiplexing of speech with data in telephone channels 
and, in certain instances, buffering of data generated by the coding of 
moving images in the Picturephone® system. The switched system con- 
sists of a uniform source, buffer, switch and channel. The source feeds 
data to the buffer at a uniform rate. The buffer's access to the channel is 
controlled by the switch; if the switch is closed, the buffer empties to the 
extent of the channel's transmission rate. The on-off pattern of the switch 
is indicated by a — 1 burst process \Ej), j = 0, 1, 2, ■ ■ ■ ; if Ej = 0, 
the switch is closed for the duration [j , j + 1 ) . The burst phenomenon is 
introduced to account for two different processes responsible for the event 
Ej = 0. There are relatively long periods during which Ej = uniformly, 
and the activity separated by such periods is defined to be a burst. During 
a burst, Ej = only infrequently. The duration of a burst is an inde- 
pendently distributed random variable with a geometric or weighted sum 
of geometric distributions. The inter-burst periods are assumed to be 
sufficiently long for the buffer to empty at some point during these periods 
of inactivity. During a burst \Ej\ is a Bernoulli sequence of independent 
random variables. 

Exact expressions for a variety of performance functional^ related to 
the system described above are obtained, together with qualitative results. 
Recursive formulas are obtained for the following : (*) steady-state distri- 
bution of buffer content for a finite buffer of size N ; (ii) mean time for 
first passage across a level N ; {Hi) the probability of overflow, for a given 
level N , during a burst; (iv) mean time for first passage across a level N 
during a burst. The recursion in each case is with respect to N. The 
asymptotic behavior of the main recursions is determined. 

* The sequence of names was determined by coin tossing. 
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1. INTRODUCTION 

A convenient framework for an unified analysis of a variety of 
digital communication systems involving buffering-some are discussed 
later-is provided by the system in Fig. 1. The source emits data uni- 
formly at the rate of one symbol per unit time. The transmission rate 
of the channel is (A; + 1) symbols per unit time where A; is some positive 
integer. The buffer has access to the channel only when the switch 
is closed. The switch is controlled by a burst process [Ej], j = 0, 1, 

2, ■ ■ ■ . Ej, for every j, is either or 1. If Ej = 0, the switch is closed 
for the duration [_j, j + 1) ; otherwise the switch is open. The burst 
process is introduced to account for cases where two basically 
different types of phenomena are responsible for the event Ej = 0. 
There are relatively long periods during which Ej = uniformly; 
the activity separated by such periods is defined to be a burst. On 
the other hand, during a burst, Ej = only infrequently. The dura- 
tion or length of a burst is a random variable. It is assumed that 
the burst length is independently distributed with a geometric or 
a weighted sum of geometric distributions. The interburst periods 
are assumed to be sufficiently long for the buffer to empty during 
these periods. The statistical assumption made in the paper about 
the controlling sequence [Ej] within a burst is that it is a Bernoulli 
sequence of independent random variables and Pr {Ej = 1 } = ir where 
< 7r < 1. In a companion paper, the case where [Ej] is first-order 
Markov will be considered. 

Important aspects of various digital communication systems are 
subsumed within the framework of the system described above. 
Diverse schemes for multiplexing data with speech on telephone 
channels 1 - 2 are representative of one class of such systems. A summary 
of the main features of the system which has been described in some 
detail in Ref. 1 follows. The central idea is to utilize the telephone 
channel during the periods of silence in speech which amount to as 
much as half of the total conversation period to transmit digital data. 
The speaker needs to have priority for the use of the channel since 
otherwise the quality of speech is impaired. Ej = 0(1) corresponds to 
the decision that silence (speech) exists during the interval [j, j + 1) 
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Fig, 1— Switched communication system. 
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so that only after it is decided that the speaker is silent does the 
buffer have access to the channel. An excellent example of the burst 
phenomenon may be found in speech monologues. Due to the presence 
of phrases in speech, two types of silences, interphrase and intraphrase 
silences, exist; 3,4 the former type consists of silences no less than 250 
ms long and this is substantially longer than the mean duration of 
(uninterrupted) intraphrase silence. 

There exists another class of digital communication systems com- 
posed of systems with only one source which transmits at a nonuniform 
rate. Most of the time the source rate is less than, say, r bits per unit 
time and r is less than the channel rate r. Occasionally, for short 
periods of time, the source rate spurts to a level r x which exceeds r. 
During such periods buffering becomes necessary. These occasional 
bursts of overloading of the channel are indicated by the [Ej\ process. 
The relation to the switched communication system of Fig. 1 is clear 
if (n — r) is normalized to unity, and (r — r ) corresponds to k. 

An example of such a system for which the analysis of this paper is 
relevant arises in buffering of data generated by the coding of moving 
images in the Picturephone® system. 5 In this case, of course, r and 
n should be interpreted as average rates 6 in the two regimes, or, when 
the worst case is of interest, as the extreme rates. The results of this 
paper appear to be relevant 7 for variable rate in-frame coding, since 
during bursts of high detail, the correlation of the data rates for 
successive picture elements is not high. For frame-to-frame coding 
the first-order Markov model, to be treated in a companion paper, is 
of interest. 

Exact expressions for diverse performance functionals related to 
the system in Fig. 1 are obtained, together with qualitative results. 
As a whole they provide a rather comprehensive set of criteria for the 
design of the important parameters of the system, such as the buffer 
size and the transmission rate of the channel. A summary of the main 
contributions follows.* 

(i) A recursive formula is obtained for the steady-state distribution 
of buffer content for finite buffers. The recursion is with respect to N 
where N is the size of the buffer. 

(u) It is proved that F N , the mean time for first passage through a 
level N, is given by 

1 1 - 7T 1 

Fn = —Fn-i Fjv-fc-i H 



* N is used to denote both the buffer size [as in (i), (Hi) and (vi)"] and a level [as 
in (ii) and (*»)]. In what follows, the definition of .V should be clear from the context. 
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(Hi) gN, the probability of overflow of a buffer of size N during a 
burst, the duration of which is distributed geometrically with a param- 
eter p is given by 



/ 1 1 1 - 7T IV 

g N = I j 

\pir gN-\ T gN-k-u 



(iv) A closed expression and a recursive formula are obtained for the 
mean time for first passage through a level N during a burst; the 
recursion is with respect to N. 

(v) The asymptotic behavior of all formulas in items (i) through 
(Hi), as N becomes large, is given. The derivation is dependent on the 
following: of the roots of the polynomials associated with the recur- 
sions, either one or two roots, depending on which recursion is being 
considered, lie outside the unit circle. 

(vi) It is proved that under certain conditions on the initial probabil- 
ity distribution of the contents of the buffer, the probability of a 
buffer being full is a monotonic, nondecreasing sequence with respect 
to time ; if the buffer is initially empty, the above-mentioned conditions 
are satisfied. One of the main implications of the result is that the 
steady-state probability of the buffer being full is an upper bound 
on the probabilities of the buffer being full at any instant. Furthermore, 
a particularly simple recursion is obtained for Pn, the steady-state 
probability of a buffer of size N being full : 



/l 1 1 - TT IV 
Pn = ( ) 

\W Pn-1 TT Pn-Ic-U 



(Observe that 1/Pn is also the mean time for recurrence of the state 
corresponding to a full buffer.) 

The closed expressions obtained are for all k and N, and the recur- 
sions hold for all N s£ 2k + 1. Wherever applicable, the buffer is 
assumed to be initially empty. An important feature of the given 
formulas is that they are also given in the form of recursions. The 
advantages of recursive formulas over the alternate versions need to be 
emphasized. For a given N, typically, a closed expression for a recursion 
involves inverting a matrix of order N. For large N, the effort is 
substantial. If, in addition, it is borne in mind that a designer is 
interested in functionals associated with a range of possible buffer 
sizes, the advantages of recursive formulas of the form given in this 
paper become overwhelming. This is only to be expected since the 
recursions are obtained by taking into full account the structure of 
the matrices involved. 
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II. EQUATIONS OF PROCESSES 

Let Bj be the number of symbols in the buffer at the jth instant. 
For a finite buffer of size N, 

B J+1 = Ma,x{Bj - k, 0} if Ej = 
= Min{ Bj + 1, N\ if Ej = 1. 

Since Bj depends only on Bj-i and £/_j, the state of the Markov chain 
of interest at the jth instant, Sj, is determined by the value of Bj 
where Bj £ {0, 1, 2, • • • , N}. Let P m (n) denote the probability of the 
state S m = n. Then 

P.(0) = (i - tt) E p B _ 1 (i) (i) 

P m {i) = rPn-iii - 1) + (1 - TT)P m -i(i + k) 

i= 1,2, -'-,N -k (2) 

P.W = *P m -x(i - 1) 

i = N - k + 1,JV - k + 2, >>-,N - 1 (3) 

P m (N) = 7r[P m _!(iV - 1) + P m -i(Nn (4) 

It is well known from the theory of Markov chains 8 that the limiting 
distribution of the states P{i) is obtained from (1) through (4) by 
equating P m (i) and P m -\(i) to P(i). 

2.1 Equations for Some New Probabilities 

Central to most of what follows are the probabilities Q m (i), where 

Q m (i) = Prj (S m = i) O (Bj ^N, j ^ m)\ 

and the buffer size exceeds N. For convenience, let X m denote the 
event Sj G {0, 1, 2, • • •, N\ for all j, ^ j ^ m, so that 

Q m (i) = Vr{(S m = i)f\X m \. (5) 

The equation governing the transitions of {Q,} is derived. It is shown 
that there exists a matrix A which relates [Qi] to (Qt-ij, i.e., 

QiU) = E A*Qi-i{i) (6) 

2=0 

or, in matrix notation, Q, = AQi-\. 



11 



THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1973 



In (5) iE {0,1, ■•-,#}, so that 

Q m {i) = Pr{(S m = oni^i 



Hence, 






Q-(0 - L Pr{(5„ = i) | (S_i = j) O X ra _!}Pr{ (5 m _i - j) H X m _i} 

j = 

= LPr|(S m = i) | (S m _i = j) n X«-i}Q«-i( j) 

i-o 

(1 - r) L Q«-i(i) if t - (7a) 

*-Q„_i(i -°1) + (1 - *-)Q*-i(« + fc) 

if i = 1,2, ...,JV - ft (7b) 

»Q«_i(t - 1) 

if i = tf - A; + 1, N - k + 2, • • • , N. (7c) 

(7) defines the (N + 1) by (N + 1) matrix A. Sometimes when the 
need arises, the (N + 1) by (N + 1) matrix A associated with a given 
N will be specified by A(N). 



A = 



1 2 A; + 1 k + 2 

(1 - r) (1 - tt)--.(1 - w) 

T (1 - *■) 



N + 1 



■ (1 - r) 



•IT 



JV - k + 1 



AT 

JV + 1 



(8)« 



It will be observed that the only difference between (7) and the 
transition equations (1) through (4) for a finite buffer of size N, is 
that eq. (4) is modified since a transition from state N to state N is 
not possible in the present context. For the same reason, the matrix A 
is not a Markov matrix since the sum of the elements of the last, 
i.e., (N + l)th column is (1 — x), the remaining columns sum to 
unity as is the case for all columns of Markov matrices. 

* The dots indicate continuation of the values of the adjoining elements; remaining 
elements are assumed to be 0. 
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If the transition matrix of the basic Markov process, i.e., the matrix 
defined by eqs. (1) through (4), is irreducible, then (I — pA) is 
nonsingular for \p\ ^ 1. The proof follows from a well-known result 
in matrix theory 9 which in this case states that if all the columns of 
(I — pA) are weakly column-sum dominant and at least one column 
of (I — pA) is strongly column-sum dominant, then the matrix is 
nonsingular. 

III. STEADY-STATE PROBABILITIES FOR FINITE BUFFERS 

In this section, a formula is given for recursively generating the 
steady-state probabilities P{i) where the recursion is, with respect to 
N, the size of the buffer. To distinguish the steady-state probabilities 
for different buffer sizes, the symbol P N (i) is introduced to denote 
P(i) for a buffer of size N. 

If N ^ k + 1, as is almost always the case, an equation of the type 
given in (2), namely, 

1 1 - 7T 

P N (i - 1) - -P"(i) + P N (i + k) = (9) 

■K T 

occurs at least once and since N ^> k usually, the main body of equa- 
tions defining the steady-state probabilities is of that form. It is 
proved in Ref. 1 what may reasonably be expected, namely, every 
solution of the homogenous set of equations that define the steady- 
state probabilities is of the form 

fc+i „_ . 

P N U)= E6.M. ' 3=0,1, ■••,# (10) 

where m are the simple roots of the polynomial 

1 1 - 7T 

n k+1 - -p. k + (11) 

IT IT 

If the polynomial has multiple roots the obvious modifications must 
be made. [Note: Since < ir < 1, the polynomial in (11) has distinct 
roots whenever -k ^ k/(k + 1) ; when t = k/(k + 1), the only repeated 
root is 1.] 

The complete recursive formula for P N (j) is obtained in two parts. 
First, a recursive formula for a set of solutions Qn(J) to the steady-state 
equations is obtained and, second, a recursive formula for the 
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normalizing constant I,n is obtained. Finally, 

P"U) = -=-i*U) j - o, l, • • •, N. (12) 

3.1 Recursions for {qN(j)\ 

Let 

q N (N) = 1 (13) 

and suppose {qnU)} satisfies the steady-state equations of a finite 
buffer of size N. Hence, q N (j) has the form given in (10).* For fixed N 
and i — 1, 2, • • • , k + 1, let 

di = L *ms ■ 

7-1 

The transformation [a,-] — * (d*} is invertible since the Vandermonde 
matrix is nonsingular. Now, 

tt? w-(iir-*fi) 
y-i 
- tfir(tf -t + 1) i = 1,2, ...,Jb + l. (14) 

Also, from the steady-state equations themselves, 

di = <?*(#) - 1 (15) 

1 - 7T 

di = q N (N - i + 1) = — — » = 2, 3, • • •, k + 1. (16) 

7T* l 

Hence, significantly, {d,} is independent of N from which it follows 
that { a, } is also independent of N. 

... * +1 ifHH 



ft+i 

= E «. 

i-l 



1 (N+k)-j 1 — 7T at_; 

-Mi M< 

7T 7T 



from (11) 



1 1 - 7T 

= -ff*+*(i) QifU) J = 0, 1, ■ • •, tf. (17) 

7T 7T 

The formula for [qN+k+i(j)\ is complete if (15) and (16) are appended, 

* To distinguish between \P N (j) | and [qirtf)}, denote the coefficients in the form 
for the latter by (a,}, i.e., bi = (l/2Ar)a,-. 
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i.e., 

1 - 7T 

q N+k+ x(N + i) = -— - »-l,2,..- f * (16) 

q N+k+1 (N + A; + 1) = 1. (15) 

3.2 Recursion for the Normalizing Constant 

Let 

2at = £ fljr(j). (18) 

y-o 



Now 



JV+fc+1 k /]\» 

£ g w+H i(i) = 1 + (l -») E (-) 

1 
fc* 



(19) 



7T 

Summing both sides of (17), 

1 lr 1 "I 1 - 7T 

Sat+a+i 1 = ~\ ^N+k — 2at 

7T* TL a , * +1 J 7T 

i.e., 

1 1 — 7T 

ZAT+fc+l — — Sat+/c — Sjv. (20) 

7T 7T 

(20) is the recursion for the normalizing constant. The derivation of the 
recursive formula for \P N (j)\ is now complete. 

Observe that in the course of the above analysis, a simple recursive 
formula for the rather important steady-state probability of the buffer 
being full, i.e., P N (N), has been obtained. 

qs(N) 1 
P*(N) = i^^ = — (21) 

2j\r Sat 

and Siv, of course, is obtained from (20). 

IV. MEAN FIRST PASSAGE TIME 

Suppose N is a fixed positive integer and the buffer capacity is 
greater than N. A functional that provides substantial insight into 
the problem of designing a buffer for which the probability of overflow 
is small is F N , the mean time required for the buffer contents to first 
exceed N. It is particularly useful in the context of burst processes 
where only incomplete data are available concerning the burst length 
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distribution-provided that the length of bursts is bounded, a simple 
comparison of the bound with F N provides an useful guide. In this 
section a recursive formula for F N , the recursion being with respect to 
N, is obtained. To correspond with the practical situation, the buffer 
is initially assumed to be empty ; the same recursive formula holds for 
the other interesting initial condition, namely, the buffer initially 
contains an unit symbol. 

X m is the event that S,- E {1,2, ■ ■ -, N] for all j, g j | m. 

Oi = Pr{ overflow occurs for the first time at i] (22) 

= Pr{ (Ei-! = 1) O (S^ = N)C\ Xi-t) 

= 7r Pr{ (Si-i = N C\ Xi_i)}, from the independence of [E{\ 

= irQi-tW) (23) 

where {Qi} is as defined in eq. (5). It has been shown in Section 2.1 that 

Qi = AQi-x. (6) 

Hence, 

Oi = vQi-xiN) 

-r(0, •••,0,1)Q,_ 1 
= tt(0, .--,0, lM^Qo 

= ireU^Qo (24) 

where e, denotes the vector* with a single element equal to unity at the 
tth location and all remaining elements 0; r = N + 1. Finally, 

Fn = Mean time for first passage through level N 

CO 

- L iOi (25) 

t=i 

CO 



= *rerf L*A<-MQo 



Let 



= ire T (I - A)~ l (I - A)- 1 Q - (26) 



A 



x l = e l (I - A)' 1 



* The superscript t denotes the transpose of a vector. 
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so that 

x'(I - A) = e r - (27) 

But 

x' =-(1,1, -.-.l) (28) 

is a solution of (27) since the elements of the last, i.e., (N + l)th 
column of A sum to (1 — ir) and the remaining columns sum to 1. 
Moreover, (28) is the unique solution of (27) since (I — A) is non- 
singular. Hence, 

F N = (1,1 •••,l)(/-^)- 1 Qo 

= l'(7 - A)~ l Qo (29) 

where 1 denotes the vector with all elements equal to unity. In the 
following section, the above formula with Q = ei is analyzed further 
to yield a recursive formula. 

4.1 Recursive Formula for Fn 

Henceforth, it is necessary to be specific about the dimensions of 
A -the matrix A associated with a given N is denoted by A(N). 
The buffer is assumed to be initially empty, i.e., So = or, equiv- 
alently, Q = e\. 

Since* \I — A(N) \[_I — A(A0] -1 ei is the vector of (signed) 
cof actors of the 1st row of [/ — A (A/)] 

1 7 - A(N)\l'lI - A(N)J- l ei = \D(N)\ (30) 

where D(N) is the (N + 1) by(]V+l) matrix obtained from A(N) by 
replacing all elements of the first row of A (N) by unity. Then, from 
(29), 

\D(N)\ 

Fn = -— L - i —■ (31) 

\I-A(N)\ 

Adding rows 2, 3, ■ • -, (N + 1) of [7 - A(N)2 to the first row, 
it can be verified that 

\I - A(N)\ = t n+1 . (32) 

In Appendix A it is shown that 

\D(N)\ = \D(N - 1)| - (1 - r)ir k \D(N - k- 1)| +t* (33) 

* | A' | denotes the determinant of the matrix A". 
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Hence, 

\D(N)\ __ 1 \D(N-1)\ (1 -it)t* \D(N-k- 1)| 1 

i.e., 

1 1 ~ 7T 1 

F* - -Fjf.! JV_*_i + - (35) 

ir ir ir 

The above relation is the desired recursive formula for the mean first 
passage time. It was obtained under the assumption that the buffer 
is initially empty. An alternative assumption about the initial distribu- 
tion, which is also of interest, is that the buffer contains an unit symbol, 
i.e., aS = 1. It may be shown that even for this case the mean first 
passage time satisfies the formula (35) though, of course, the initial 
conditions to the recursion in the formula are different. 

V. PROBABILITY OF NO OVERFLOW IN A BURST 

The results of this section are useful when information concerning 
the length of bursts is available. It is assumed that the distribution of 
burst length may be expressed as a weighted sum of geometric distribu- 
tions. Given below are formulas which yield the probability that the 
contents of the buffer during bursts do not exceed N, a given positive 
integer. 

At this stage, assume that the distribution of burst length is geomet- 
ric; the generalization to distributions that are weighted sums of 
geometric distributions will be taken up later. If the burst length is 
denoted by I, then 

p r {j = i) = (i - p ) p i-i »-l,2, .... (36) 

for some p, < p < 1. Let Gn = Pr {buffer contents do not exceed N 
during a burst}. The usual decomposition into mutually exclusive 
events yields 

G N = £ Pr{Sy G (0, 1, • • -, N), j = 0, 1, • ■ •, m; 

and burst length = m\ 

= L Pr{X m nZ = m) 

= L Pr{X m }Pr{Z = m). (37) 

m>l 
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The last relation holds since, Pr{X m \l = m\ = Pr{X m }. Now, 
Pr{X m } = lPr|S m = ini n | 

N 

= £ Q.(0 from (5), 

i=0 

- VQ m 

= VA m Q (38) 

where A is the (JV + 1) by (N + 1) transition matrix denned in 
Section 2.1 and Q is the vector given by the initial distribution-it 
may be assumed that S G (0,1, • • •, N). Hence 

Gx = L l'A"Qo(l - p)p m ~ l 



P 
1 -p 



L (M)' 

m>l 



Q< 



!'{(/- pA)-i-/}Qo 



p 
(1 -p) 



|l'(J-pA)-»Qo-l}- (39) 



In the sequel, a recursive formula for Gn is developed for the case 
where <So = or, equivalently, Qo = &\- 

5.1 Recursive Formula for G N 

The matrix A associated with a given N is denoted by A(N). 
\I — pA(N)\ {I — pA(N)}~ l e\ is the vector of (signed) cof actors of 
the first row of {I- P A(N)}. Therefore, \I-pA(N)\V[l-pA(n))- l ei 
is the determinant of the matrix B(N) obtained by replacing every 
element of the first row of \I — pA (N) | by unity. 

\B(N)\ 
I'll - pA{N)\-^ = ' ' — • (40) 

\I - pA(N)\ 

Let the (signed) cofactor of the element }/ — pA(N)]u be denoted by 
C", i = 1, 2, • ■ -, N + 1. From the definition of B(N), 

N+l 

\B(N)\ = £ C". (41) 
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The elements of the (N 4- l)th column of {I — pA(N)} sum to 
{1 — p(l — 7r)| and the elements of the remaining columns sum to 
(1 — p). Hence, by adding the rows 2, 3, ■ • • , N + 1 to row 1, it follows 
that 

\I - P A(N)\ = (1 -p)EC lf + {1 -p(1 -r)}&>** 



i.e., 



N+l 

= (i -p) e c«+ {l-pd-TO-a-pMc 1 -^ 
i=i 

= (1 - P )\B(N)\ + pxC 1 -^ 1 from (41) , 
I B (AT) | 1 p7T C 1 -'^ 1 



|7 — pA(iV)| 1-p 1 -p \I - P A(N) 

Recapitulating, 

1 -P 



(42) 



f, - -[1'{1 ~ piWI^i - 1] from (39) 

B(N) 



p 
1 ~P 



r gwj ii 

l\I-pA(N)\ J 



= 1 - 7T 



|/-pA(iV) 
C'.^v+i 

\I- P A(N) 



from (40) 



from (42). (43) 



The remainder of the derivation is in two parts. First, a closed form 
expression for C 1,Ar+1 is obtained. The second part is on the recursive 
formula for \I — pA(N) \ and this formula is derived in Appendix B. 



1 2 3 k + 1 k + 2 
{ (1 +\) X X... X 

M 1 X 



A/ + 1 



/ - P A{N) = 



1 



M 1 
M 1 



(44) 



N -k + l 



N 
N + l 
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where — p(l — ir) = X and — p?r = p. C 1,lf+1 , the (signed) cof actor to 
{I — pA(N)} i,n+i, is the signed determinant of an upper triangular 
matrix ; 

= p V. (45) 

In Appendix B it is shown that if xn, a scalar, is used to denote 
\I — pA (N) | , then the following recursive formula holds : 

x N = x N -i - p fc+1 7r*(l - 7r)XAT_ A _i. (46) 

Hence, 

*n 1 / xpf-i \ 1 — t / x N -k-\ \ 

Let 

a x * \I - P A(N)\ 

y« - ^ - ^ (^ 

so that, 

1 1 - 7T 

Vn = — Vn-\ ftr-ft-i. (49) 

7Tp 7T 



.N+lpN 



From (43) and (45), 

Gjr - 1 - ■ 

|7- P A(iV)| 

i.e., 

Gfjv = 1 — — • (50) 

Vn 

(49) and (50) together provide the desired recursive formula for the 
probability that the contents of the buffer does not exceed a given level 
N during bursts if the buffer is initially empty and the distribution of 
burst lengths is geometric. 

Suppose the distribution of burst lengths is the weighted sum of 

geometric distributions; i.e., 

j 
Pr{ burst length = i\ = £ a,-(l - Pi )(p y )'- 1 . (51) 

It may then be shown that G N = Y, J j=i<xjG].n where Gj. N is obtained 
from (50) and (49) with p replaced by p, in the latter equation and 
Gj,n identified with Gn- 
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VI. MEAN TIME FOR FIRST PASSAGE IN A BURST 

In Section IV certain formulas for the mean time for first passage 
across prespecified levels are given. In burst processes where data 
regarding the length of bursts is available, a more meaningful functional 
is one in which a level is defined to be crossed only if this event occurs 
during a burst. Bursts, then, may be visualized as a period of observa- 
tion of the buffer. First passage across N, a positive integer, is defined 
to occur at i if 

[S s ^ N, j - 0, 1, 2, • -, i - 1 and S t = N + 1 and, 

burst length ^i). 

Let Ri denote the probability of this event. The functional of interest 
is H N = Li°=i iRi- The burst length distribution is assumed to be 
geometric ; generalization to larger classes of distributions may be un- 
dertaken as indicated in the preceding section. Hence, if I denotes 
burst length, 

p r {Z = i\ = (l - p) p «-i i = 1, 2, • • • (52) 

for some p, < p < 1. 

In the notation of Section 2.1, 

r { = Pr{Si-i = n n x,_i n £,-_! = mi^i\ 

= Pr{#,_i = l}Pr{Si_i = NnX^l^ i}Pr{l^ i\ 
= 7rPr{S,_i = NC\ X,-_!}PrU ^ i\ 
= TQ i -i(N) P i ~ 1 

= 7re r (pAy- 1 Qo. (53) 

A is, of course, the (N + 1) by (N + 1) matrix defined in Section 2.1 
and Qo is the initial condition vector. 

H N = Te T (j:i(pAy- 1 )Q 

i.e., 

H N = ire T {I - pA)~HI ~ pA)-'Qo. (54) 

The above concludes the derivation of the closed formula for H^-the 
rest of the section is concerned with recursive versions of the formula 
for the case where the buffer is initially empty, i.e., Qo = ei. Once again, 
it is necessary to revert to the use of the symbol A (N) to denote the 
matrix A associated with N. 



For fixed N, 
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z(p) = e r E P «A'(iV)e, 



Hence, 



= perZ { P A(N)\'ei 

= peAI - pA(N)}-i ei . (55) 

z{ P ) = -z(p) = e T £ (» + Dp'A'W* 

dp »£0 



= erZ {{pAiN}}^! 


1 
= -#*• 

7T 


Returning to z(p) and (55), observe that 


p Ci.Ar+i 

-fnl — 




* \P 1 

\I- P A(N)\ 


C 1,N+1 = pNlr N, 


Hence, 

~(n) — 




<■• \P) 

\I- P A(N)\ 


and, from (56), 


d f p N + X TT N+l 

H — 




d P \\I - P A(N)\ 
Let 

1 A ptf+lT^ 1 





t* |/-aAW 



(56) 

(57) 
(45) 

(58) 



(59) 



Since v N = {l/pir)y N where yN has been defined previously in (48) and 
the recursion in (49) for y N is linear, vn satisfies the same recursion. 
Hence, with un = (d/dp)v N (p), the following formula is obtained: 

H N = (60) 

2 

1>N 
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and, 

Vn = Vff-i VN-k-l (61) 

pir t 

1 1 - T 1 

Uff = Utf-l u N -k-i — Vn-1 . 

pir t p'tt (OJJ 

VII. ASYMPTOTICS OF RECURSIONS 

The main recursions occurring in the paper are of the following 
forms : 

x N = - Xat-i XN-k-i (63) 

IT T 

1 1 ~ * 1 

VN = - Vn-i • VN-k-i + - (64) 

7T 7T T 

z N = — zn-i ZN-k-i where < p < 1 . (65) 

pir tv 

Equation (63) occurs in the formula for the (unnormalized) steady- 
state probabilities and in the formula for the normalization constant ; 
(64) occurs in the formula for the mean first passage time ; (65) occurs in 
the formula for the probability of no overflow during bursts. The 
fundamental solutions of these recursions are obtained from the roots 
of the following polynomials. 

F{p) £ M *+' - - n k + ^^- (66) 

TV 7T 

G(n) = y.«+i M * + -^- (67) 

pir it 

Equation (66) is associated with (63) and (64); (67) with (65). The 
two results given below enumerate and estimate the roots of F(p.) and 
G(p) outside the unit circle. 

Lemma l' : Except for one positive real root 1/6, and 1, all other roots of 
F(n) lie inside the unit circle |/i| ^ 1. The root 1/6 lies outside the unit 
circle if and only if k > ?r/(l — ir). 

Lemma 1 is a specialization of a result proved in Ref. 1. Bounds on 
6 are also given there. 

Lemma 2: G(/x) has k roots inside the unit circle |/t| ^ I, no roots in the 
annular ring 1 ^ |ju| ^ 1/p, and one real, positive root outside the 
circle \p,\ ^ 1/p. 
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Proof: 

1 - 7T 

G(0) = > 

7T 

G(l) - - (1 - 1/p) < 

7T 

1 - ir r 11 

w L p a+1 J 

Since G(0) > and G(l) < 0, there exists a real positive root of 
G(ji), r, where r < 1. Since G(l/p) < and G((i) -» °o as /i — > oo, 
there exists a real positive root of G(p), R, where i? > 1/p. The 
following theorem which is stated without proof may now be applied. 

Pellet's Theorem: 10 Given the ■polynomial 

f(z) = a + oiz + • • ■ + a p zp -f • • • -f a„z n , a p ^ 0. 

7/ the polynomial 

F p (z) = \a \ + |a x |2+ ••• + |a p _i|z p - 1 

- \a p \z p + |o p+ i|2P +1 + ••• \a n \z n 

has two positive zeros r and R, r < R, then f(z) has exactly p zeros in or 
on the circle \z\ < r and no zeros in the annular ring r < \z\ < R. 

Identifying p with k, n with k + 1 and f(z) with G(p.) the rest of 
the proof follows. 

The reader may now verify that, for large N, 



~C if lb ' * 








1 - 7T 




S C, + C 2 N if k = tt/1 


— 7T 


9* Ci + C 2 ( - ] if k > - 


IT 



1 - 



2/^C, + JVC 2 if k < 



1 - 7T 

S* C x + JVC, + N 2 C 3 if fc = tt/1 - 7t 

7T 



1\tf 



» C, + #C 2 + C,( - if A; > 



1 - T 
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where, R and 1/0 are roots previously denned and the C's are constants. 
The constants may be obtained by fairly straightforward computations. 
The qualitative difference between the forms of the expressions 
corresponding to k < ir/(l — *■) and k > ir/(l - t) are noteworthy. 
This is not unexpected, since it may be recalled that in Ref. 1 it was 
proved in a more general context that the Markov chain associated 
with the infinite buffer is positive recurrent if and only if k > v/(1—t). 

VIII. A MONOTONICITY PROPERTY OF THE PROBABILITY OF A FINITE 
BUFFER BEING FULL 

The steady-state probability of a buffer being full, i.e., P(N) where N 
is the size of the buffer [see Section III and, in particular, eq. (21)] 
may be expected to be an important factor in the practical design of 
buffers. This is so not only because of the immediate implications of 
the definition but also because l/P(N) is the average recurrence time 
of state N. However, this approach would appear to overlook the 
possibility that the probability of the buffer being full in the transient, 
i.e., in the approach to steady state, is seriously underestimated by 
P(N). Such an event is not easy to rule out because, after all, P(N) is 
an element of only one (normalized) eigenvector of the transition 
matrix while all the modes or eigenvectors and eigenvalues of the 
matrix contribute to yield P m (N) when m is finite. However, one of the 
implications of the result in this section is that, under certain conditions 
on the initial probability distribution of the contents of the buffer, 
P(N) is indeed an upper bound on P m (N), i.e., P m (N) ^ P(N), 
m = 0, 1, •• • ; furthermore, the important case of the buffer being 
initially empty satisfies the conditions just mentioned. 

For a buffer of size N, the result states the following. Suppose at the 
with instant the state probabilities satisfy the inequalities : 

irP m (i) - (1 - tt) £ P m (j) ^ i - 0, 1, • • •, N - k (68) 

JV 

rP m (i) - (1 - tt) £ P m (j) ^ 

i = N - k + 1, N - k + 2, • • •, N - 1. (69) 

Then (a) P m (N) ^ P m +i(N), and, as shown below, (b) the inequalities 
in (68) and (69) are satisfied with P m (l) replaced by P m +i(l) for I = 0, 
1, 2, ■■■ , N. Therefore, if (68) and (69) hold, Pi(N) g P i+l (N) for 
all i f i ^ m ; i.e., the probability of the buffer being full is a monotonic, 
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non-decreasing sequence, (a) may be trivially verified. The proof of 
(b) is as follows. 

(i) i = 0. 

irP m+1 (i) - (1 - w){P m+ i(i + 1) + P m+1 (i + 2)- • • + P m+l (i + ft)} 

= 7T(1 - 7T)[P TO (0)+ P m (l) + • • • + P m (ft)] - (1 - 7r)7T 

X [P B (0) + P»(l) + • • ■ P m (ft - 1)] - (1 - 7r) 2 [P m (fc + 1) 

+ P m (k + 2) + • • • + P m (2fc)] 
= (1 - 7r)[7rP m (fc) - (1 - w){P m (k + 1) 

+ P m (fc + 2) ••• + P m (2ft)|] ^ 

(w) 1 g * ^ iV - 2ft. 

irP«+i(t)- (l-ir){P m+ i(i+l)+P m+1 (i + 2) ••• + P m+1 (i + ft)} 

= 7r(l - 7T)[P m (0) + P m (l) + • • • + P m (ft)] - (1 - 7T)7r[P ro (0) 
+ P.(l) + • ' • + P m (fc - 1)] ~ (1 - 7r) 2 [P m (ft + 1) 

+ P m (k + 2) + ■ ■ • + P m (2ft)] 
= (1 - 7r)[irP m (ft) - (1 - 7r)|P m (ft + 1) 

+ P m (ft + 2) + ••• + P„(2fc)}] ^ 

(Hi) N-2k + lgi£N-k-l. 

xP m+1 (i) - (1 -7r){P m+1 (;+l) 

+ P m+i {i + 2) + ■ • • + P m+ i(i + fc)j 

= 7T[P m (z - 1) + (1 - 7T)P m (i + ft)] - (1 - t) 

X 7r[P m (z) + P m (i + 1) + • ■ • + P m (i + ft - 1)] 
- (1 - 7T) 2 [P m (i + ft + 1) + P m (i + ft + 2) + • • • + P m {N)~\ 
= r[rP,(* - 1) - (1 - *-){P B (i) + P n (i + 1) 

+ ■ • • + P„(i + ft - 1))] + (1 - 7r)[>P m (t + ft) - (1 - tt) 
X |P.(i + fc+ 1) +P*(i + ft + 2) + ... +P m (A0J] ^ 

(iv) i = N — ft. 

7rP m+1 (t) - (1 -7r){P m+1 (z + l) 

+ P m+1 (i + 2) + • • • + P m+1 (t + ft)} 
= ir[>P w (i- - 1) + (1 - 7r)P m (i + ft)] - (1 - 7r)7r[P m (z') 

+ P m (i + 1) + • • • + P m (i + ft - 1)] - (1 - T)irP m (N) 
= x[*P.(*- 1) - (1 - w){P m (l) 

+ P m (i + 1) + • • • + P m (i + ft - 1)}] ^ 
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(v) N - k + 1 £ i ^ N - 1. 

rP-nffl " (1 - 7T)fP m+1 (i + 1) + P m+ l(i + 2) + - ■ 

- 7r2P m (i - 1) - (1 - »>[P.(»") + P m (i + 1) + ■ 
= 7r[TP m (i- l)-(l-*)\P m (i)+P m (i+l) 



+ P m+ i(iV)i 
■ + P.OT] 



+ ••• + P-(tf)}] ^0 

(b) is proved. 

Observe that P«(0) = 1, P«(i) = 0, i = 1, 2, ■■■ , N satisfies the 
inequalities (68) and (69). However, the other initial distribution of 
interest, namely, P m (0) = 0, P„(l) = 1, P„(i) = 0, i = 2, 3, ■ • ■ , N 
does not satisfy the inequalities. Also, it may be verified that for the 
latter set of initial conditions, the monotonicity property does not hold. 

It is interesting to note that if (68) and (69) hold, then P m (0) 
^ P m+ i(0), so that together with (b), P,(0) ^ P,-+i(0) for alii, i ^ m, 
i.e., the probability of the buffer being empty is a monotonic, non- 
increasing sequence. 



APPENDIX A 

Recursive Formula for \D(N)\ 
2 k + 1 

Or - 1) 



1 

r i i . . 

-7T 1 



D(N) = 



"1 

-7T 1 



N N + l 
, 1 1 



•1 

— IT 



Or- 1) 




Expanding \D(N) | along the (N + l)st row yields 

\D(N)\ = \D(N - 1)1 +tt|X| 
where 



N -k + 1 
JV - fc + 2 



N 

N + 1 



(70) 
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X± 



12 k + 1 

11.... 

-7T 1 (7T - 1) 



N - 1 
. 1 



N 
1 



7T 1 

-7T 1 

-7T 1 



.(tt-1) 



(tt-1) 



7T . 1 



N - k- 1 

N - k 

N - k + 1 



.V 



Expanding | X | along last column : 

|X| = (-D^+^-tt)^- 1 + (tt - l)(_l)AT+(JV-*+D|y| (71) 

where 



Y^ 



1 2 

1 1 . 

-TT 1 



ft + 1 



N - 1 



- 1 



1 ~l 1 

2 



1 

-TT 1 



"TT - 1 




-TT " 1 

-TT TV - 1 



#- fc - 1 

N -k 



Expanding | Y\ along the last (k — 1) row. 

\Y\ = (-tt)*- 1 !/)^ - Jb- 1)1. 



(72) 
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Combining (70) through (72) : 

\ D (N)\ = \D(N - 1)| - (1 - *)ir>\D(N - k - 1)|. (73) 



APPENDIX B 

Recursive Formula for \I — pA(N)\ 

I — pA{N) is given in (44). Let x n denote \I — P A(N)\. 
X = — p(l — t) and n = —pir. 



Also, 



(t) Expand \I — P A(N)\ along the last, i.e., (N + l)th row, of 

x N = Xn-i — n\X\ (74) 

where 

12 A- + 1 N 



fl + X X . 
M 1 



x£ 



x o 

X 



X 

X 



/z 1 

M OJ AT 



AT-fc- 1 

AT- ft 
AT - fc + 1 



(it) Expand \X\ along the last, i.e., Nth column of X. 

\X\ = (-l)^-*+ 1)+Ar X|F| 

-(-i) w x|r| 

where 



Y^ 



1 2 

1 +x x 

M 1 



fc + 1 

. x 

X 



i o 
n l o 

li 



X 



1 

• n' 1 



(75) 



iV — Jb — 1 
AT-lb 

AT - k + 1 



AT- 1 
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(in) Expand | Y | along the last k — 1 rows. 

\Y\ - |i*-^jr_*_i. 
Hence, 

a: at = xat-i — /* | -3T | 

-*jr-i- *(-l)*-*X|F| 

= Xn-\ - fi(-l) k - 1 \li k - 1 X N - k -i 

= x N -i - p* +1 7r fc (l - r)xir-t-i. 



(76) 

from (74) 
from (75) 
from (76) 

(77) 
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